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Abstract 

In this paper, we present a modular strategy which describes key properties of the absolute 
primary decomposition of an equidimensional polynomial ideal defined by polynomials with ra- 
tional coefficients. The algorithm we design is based on the classical technique of elimination 
of variables and colon ideals and uses a tricky choice of prime integers to work with. Thanks 
to this technique, we can obtain the number of absolute irreducible components, their degree, 
multiplicity and also the affine Hilbert function of the reduced components (namely, their initial 
ideal w.r.t. a degree-compatible term ordering). 

Introduction 

In this paper we design an algorithm whose aim is quite simple to state: 

Given a set of polynomial rational equations which define an equidimensional algebraic set of C", 
we would like to get as many information as possible on the irreducible components of this algebraic 
set. 

We can rephrase the problem in algebraic language: given an ideal (with suitable hypothesis on 
its dimension) in the polynomial ring over C[X] defined by rational generators, find all the possible 
information about its primary components. 

The problem is really simple to state and many authors looked for efficient strategies to get the 
irreducible decomposition of an algebraic set: one can see for instance 1 8 1 and the references therein 
to have an overlook of the different techniques. In many Computer Algebra Systems (CAS for short) 
you can find routines computing the primary decomposition of an ideal: the underlying algorithm is 
often the one described in [13|. Nevertheless, the problem is really challenging since the existing 
algorithms and implementations often focus on particular cases, e.g 0-dimensional ideals (see for 
instance 1. 10 J ): for more general situations, also the best implemented algorithms (for instance, the 
ones in lUl or |[T3l ) may have unsatisfying time of execution and there may be problems of memory 
allocation. In fact, the computations required to a personal computer to find a primary decomposition 
are often quite heavy. 

Our aim is to design an algorithm concerning the decomposition of an ideal which can give an output 
in a reasonable time and with a limited use of memory. 

The main computational tool that we use are modular computations, taking this technique from 
the absolute factorization algorithm for bivariate polynomials presented in 131. 
Recent papers about decomposition of algebraic sets (see for instance |[T2l . |[2T1 ) focus on getting 
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information about the irreducible components from a generic section with a linear space, namely they 
bring back the problem to the study of a 0-dimensional ideal. 

In this paper we will bring back the problem of computing a primary decomposition to the problem 
of computing an absolute factorization; this technique is in some sense "classical" ( ifTBl ). but not very 
exploited because not efficient from the computational point of view; a powerful improvement of this 
technique is in [4], where the authors avoid the use of generic projections in order to compute the 
equidimensional isoradical decomposition of an ideal, using as coefficient ring Q or a finite field of 
positive characteristic. Our approach is instead to use generic projections (by a generic change of vari- 
ables and projections on coordinate linear spaces) and exploit modular computations to move around 
the computational difficulties, preserving a lot of data concerning the absolute primary decomposition 
of the ideal. The output of our algorithm will not be the complete primary decomposition of the given 
ideal, but it will return information concerning the components, such as number, degree, multiplicity 
and, for reduced components, the affine Hilbert function. 

In Section [T] we will show that once known information about one of the primary components of 
the ideal, the same is known for other components too. We simply rephrase the definition of "con- 
jugacy" for absolute factors of a multivariate polynomial with rational coefficients (see |5 |, Lemma 
9.0.8), for primary components of an ideal generated by polynomials with rational coefficients. De- 
gree, multiplicity and affine Hilbert function are "invariant by conjugacy", so if we obtain this in- 
formation about a primary component, we actually have the same information for all the primary 
components in its "conjugacy class", avoiding to repeat computations. 

In Section |2] we show that, fixed an algebraic extension L of Q, there are infinite prime integers 
that implicitly define a homomorphism from L to Z/pZ (more precisely, an inclusion of L in Qp, 
Lemma 1231) . This means that with a careful choice of a prime p, we can reduce the coefficients of a 
polynomial in L[X] modulo p. Furthermore, infinite prime numbers preserve interesting properties of 
an ideal in L[X], namely the initial ideal with respect to some degree-compatible term ordering and, 
as a consequence of this, the affine Hilbert function. 

Summing up, we can choose a prime p which allows modular computations in L (we can choose it 
using Lemma 123] ). Only a finite number of primes p does not preserve the properties of the primary 
components we are interested in, so we can assume that we are avoiding them by taking a "generic" 
prime p. 

In Section [3] we present the exact strategy to obtain the prime components of an ideal o. This 
technique is mainly based on elimination of variables, in order to bring back the problem of primary 
decomposition to a problem of factorization. This strategy was first investigated by Crete Hermann 
in 1(161 and it is similar to the splitting techniques presented in f4l, but we present it completely for 
lack of an accessible reference on the whole strategy. Nevertheless, the technique of Section[3]is not 
efficient from a computational point of view: first of all, projections are actually computed with a 
generic change of coordinates and an elimination of variables performed by a Groebner Basis; then, 
in order to obtain the reduced primary components, we compute a colon ideal; this is performed again 
by an elimination Groebner Basis. 

In Section|4]we try to gain in computational efficiency, even if we "lose" the exactness of Section 
121 We will apply the modular results of Section |2l on the exact algorithms of Section |3l We compute 
projections, factorizations and colon ideals modulo well-chosen prime integers; we do not get the 
reduced primary components, but we obtain an algorithm (Algorithm |3l) which can compute the initial 
ideal of the reduced components of a and give information about the non-reduced components. 

Finally, in Section |5l we test our strategy on a simple example, a complete intersection ideal in 
3 variables, getting the output of Algorithm [3l in a really reasonable time. The same ideal could 
not be decomposed by other CAS in 1 hour (because of problems with memory allocation); obvi- 
ously the comparison between our strategy and implemented primary decomposition algorithms is 
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not complete, since Algorithm [3] does not return the absolute primary decomposition of the input 
ideal. However this comparison enlightens promising performances of our strategy and this can be a 
starting point for designing an efficient primary decomposition algorithm. 

Notations 

In what follows, we will work in a polynomial ring R with coefficients in a field IfC of characteristic 
0: R = K[Xi, ...,Xn]= IK[X]. We will precise, when needed, if IC = Q, Q(a) or C. 

Given an ideal o C i? we will consider its zero set in C": W = V{a) = {P € C"|/(P) = OV/ G 

a}. 

1 Affine Hilbert Function and Conjugacy 

In this section we introduce the main definitions concerning the primary decomposition and the affine 
Hilbert function of an ideal a C R. 

We show that some of the primary components of an ideal a are very "similar" to each other, in the 
sense that given a set of generators for a primary component, we can get a set of generators for another 
primary component by means of conjugacy, just like we do for the absolute factors of a polynomial 
with rational coefficients ([5|, Lemma 9.0.8). This allows us to avoid repeating the computation of 
the affine Hilbert function for the conjugate components, since it is invariant by conjugacy. 

For all the definitions and properties concerning primary decomposition, the main reference is 
im, Chapter 4. 

Definition 1.1. A proper ideal q in a ring R is primary if the following condition holds: 

xy £ q and x ^ q =^ y G ^/q. 
Every prime ideal is obviously primary. 

Proposition 1.2 (OJ, Proposition 4.1). Let c\ be a primary ideal in R. Then p = ^/q is the smallest 
prime ideal containing q; we say that q is p-primary. 

Definition 1.3. A primary decomposition of an ideal a in R is an expression of a as a finite intersection 
of primary ideals: 

r 

a=f|qi. (1.1) 

i=l 

If moreover 

1- qi 1 qj.- 

2. the prime ideals pi = ^Jql are all distinct, 

then the primary decomposition di.il ) is said to be minimal. Any primary decomposition can be 
reduced to a minimal one (see page 52). 

Since we assume to work in a polynomial ring R with coefficients in a field, a minimal primary 
decomposition always exists. 

The factorization of a multivariate polynomial and the primary decomposition of a polynomial 
ideal are very close to each other: indeed, the primary decomposition of a principal ideal corresponds 
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to computing the absolute factorization of the generator of the ideal. So we can look at the factoriza- 
tion of a multivariate polynomial as a particular case of primary decomposition. 
Thanks to this similarity, it is natural to extend the definition of degree and multiplicity of a factor to 
a primary component. We can define them through the ajfine Hilbert function ( ifTTl . Section 5.6): 

Definition 1.4. Let a be an ideal in the polynomial ring R standard graded. 

We first define {R<i), the vector space generated by all the polynomials of R of degree < i. The 
"K-vector space {a<i) is the vector subspace of{R<i) which consists of the polynomials of a of degree 
< i. Since a<j = R<i R a, we can view the vector space R<i/a<i as a vector subspace of R/a. 

The map HF^^^ ■.'L ^'L defined by 

HF^/ai^) = d[mK{{R<i)/{a<i)) for i e Z 
is called the affine Hilbert function of R/a. 

From Definition 11.41 it is natural to define the affine Hilbert series, polynomial, dimension and 
the affine regularity index of i?/a. The definitions are similar to the analogous for the homogeneous 
case; for all these definitions and their properties, we refer to lITTl . Section 5.6. 

Definition 1.5. Let a be a proper ideal in R, consider its affine Hilbert polynomial HP^^^{t) G Q[t]. 
The degree of R/a is (dim(i?/a))! • {lcoeff{HP^^^{t)). 

We will often say "dimension and degree of a", meaning the dimension and degree of i?/o. 

Finally, once defined the degree of an ideal, we can define the multiplicity of a primary compo- 
nent. Here we state the algebraic definition, which corresponds to the intuitive idea that the multipli- 
city counts "how many times" the primary component is repeated. 

Definition 1.6. [I\2^, Definition 10] Let q G R be a p-primary ideal. Then the multiplicity of q, 

mult(qj), in p is deg(p)/ deg(q). 

We will often talk about the multiplicity of a primary component, implying that it is the multi- 
plicity in its radical. 

We now briefly recall how to explicitly compute the affine Hilbert function. 

For a polynomial ring ]K[X], we will denote with T" the monoid of monomials in IfC[X] and with 

= Xl^ • • • , ,ij S N a monomial. A term ordering ^ on is degree compatible if for any 
couple of monomials X^, X'^ 

X^ ^ X-^ ^ degX^ < degX"^. 

Once fixed a term ordering < on T", for a polynomial g G K[X], we denote with LM^{g) (or 
simply LM{g) if there is no ambiguity) the maximal monomial with respect to < appearing in g with 
non-zero coefficient. 

In the following Proposition, HFjin^ is the Hilbert function for a homogeneous ideal. 

Proposition 1.7. Let < be a degree compatible term ordering on T". For every i £ 7j, we have 
^^R/a(^) = 'E]=oHFR/LAUia){j)- In particular, we have HF^^^{i) = HF'^i^j^^^^^^{i) for all 

Proof See HIl, Proposition 5.6.3. □ 
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Proposition 11.71 gives us the practical way to compute the affine Hilbert function of a: chosen 
a degree compatible term ordering ^, we can compute the initial ideal of a and then we count the 
number of elements in the vector space R/LM^{a){j) for every j < i. 

We now show with a few lemmas that, given an ideal a defined by polynomials with rational 
coefficients, the relation of the primary decomposition on Q[X] and the primary decomposition on 
C[X] is similar to the relation between a rational and an absolute factorization of a multivariate 
polynomial with rational coefficients (as shown in [5|, Lemma 9.0.8). In other words, there is a 
conjugacy relation among some of the primary components. 

Definition 1.8. Consider a ideal in C[X] defined by a set of polynomials with rational coefficients. 
Let a = f]j=i ({i G Q[X] be the rational primary decomposition of a; qj (resp. V{qi)) is a rational 
primary component of a (resp. ofV{a)). 

We then consider a primary decomposition in C[X] of each rational primary component c\i: 

Ifri = 1, we say that (resp. V{qi)) is a pure rational component of a (resp. ofV{a)). 

Consider a non-pure rational component q^ of an ideal a C i?. Let Lj be the smallest (w.r.t. 
the degree of extension on Q) normal algebraic extension of Q such that q has a set of generators 
in Lj[X]; assume that Lj = Q(ai) and we denote with /j(aj,X) a polynomial in the chosen set 
of generators of q^; indeed, we can think of such a generator as a polynomial in Q[Z, X] with Z 
evaluated in aj. Consider the Galois group of Lj over Q, Qi = Ga/(Lj/Q). 

For every a G Gi, starting from q^-^^ with ^J^^ = pf'K we can define an ideal in the following 
way 

qf ) = (/i(ai,X), . . . , fiiai,X)) ^ a(qf = {fi{a{ai),X), fi{a{ai),X)). 

Obviously, the definition of a{ql^'^) is independent from the chosen set of generators of q^^^''; q-^^ and 
a{q[^^) have the same dimension; furthermore it is straightforward that the ideal a{q\^^) is a{p\^^)- 
primary. Finally, if r,^ G Gi, t ^ a, then T(qf ^) / a(qf ^) and r(pP) / fT(pf ^). 
We now show that actually the ideals o"(q-^^) are the primary components of q^ in Q[X]. 

Lemma 1.9. Consider q j a non-pure rational component o/ a C C [X] defined by polynomials with 
rational coefficients, Lj the smallest normal algebraic extension of Q such that q^^^ has a set of 
generators in Lj[X], G = Gal{L,i/Q). The minimal primary decomposition ofqi is 

and in particular = ^Gi = [Lj : Q]- 

Proof. The ideal o-(q^^^) is cr(p^^'')-primary. So Hcree ^ primary decomposition of an 
ideal b. 

Furthermore, it is a minimal primary decomposition. Indeed, thanks to the definition of the ideals 
through the automorphism of Lj, all the associated primes (Ti{p\^^) are distinct; for what concerns 
redundant primary components, for any a ^ G, since Lj is the minimal normal algebraic extension 
containing a set of generators of q^/'\ then there is / G q^-^^ such that Yir^a ''"(/) ^'^^ '^(l!^'')- 
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Using [1], Exercises 12 and 13 of Chapter 5, if we consider the associated primes and the natural 
homomorphism Q[X] — )• Lj[X], then the set of prime ideals {pij} is the same as the set of prime 
ideals of Lj[X] whose contraction is pj. Then Q acts transitively on the set that is 

So ricreS ^ minimal primary decomposition of q^; since is a primary ideal, all its 

primary components in Lj[X] have same dimension and there are no embedded components, so its 
primary decomposition is unique and is exactly the one in (I1.2I ). □ 



Lemma 1.10. Consider a Q C[X], defined by a set of polynomials with rational coefficients, with no 
embedded components. Then the minimal primary decomposition of a is 

-=nfn-(^.'")V 

i=l \(TGgi / 

Proof. Since the ideals a ^g^^^^are primary, we just need to show that the decomposition is minimal. 
Condition [T] of Definition 1 1 . 3 1 about minimality is straightforward from Lemma [L9l 
For what concerns Condition [2l we just have to point out that if there is p associated to two 

different primary components of o in C[X], a {q^j^^^, cr £ Qi, and r [cfp^, r G Qj, then we would 
have two associated primes of the rational primary decomposition included in one other or equal. 
But this contradicts the fact that a has no embedded components and the minimality of the rational 
primary decomposition of a. □ 

Up to relabeling the automorphisms of GaZ(Lj/Q), we can rewrite (11.31 ) as 

withri = [Li:Q],qp)=aj(qf) 

Definition 1.11. The primary decomposition f |j.4D is the absolute primary decomposition of a and for 
i such that r, > 2, we say that qp^ (resp. V{(('p)) is an absolute component of a (resp. ofV(a)). 

If qi = rijLi li"''*' we say that q^-* and q^^-^ ■* are conjugate, and that {qp-*}j=i is a conjugacy 
class. Any number or property of an absolute component is invariant by conjugacy if it is the same 
for all the absolute components in the same conjugacy class. 

From now on we will focus on a particular kind of ideals, equidimensional ones. 

Definition 1.12. An ideal a (resp. an algebraic set W C C"j is equidimensional if all of its primary 
components {resp. all of its irreducible components) have the same dimension. 

Thanks to H], Corollary 4.11, if a is equidimensional, all of its primary components are uniquely 
determined and so in this case the primary decomposition is unique. 

Furthermore, there is a wide class of ideals which are equidimensional: if we consider a complete 
intersection ideal a G i? generated by n — c polynomials, it is equidimensional. This can be seen as a 
consequence of the Affine Dimension Theorem ( |[T5l . Chapter I, Proposition 7.1). 

We can finally fix our purpose. 
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Given a non-prime equidimensional ideal a C C[X], generated by polynomials witli 
rational coefficients, we write its primary decomposition as in (I1.4I ). Then for every 
rational primary component of o, we would like to find the numbers r^, deg(qj-^'')and 
mult(qf^). 

If mult(qj-^^) = 1 (the primary component is reduced, and so it is prime) then we would 
also like to compute the affine Hilbert function of i?/q-^\ 

Remark 1.13. Thanks to Lemma {1.10\ all the information concerning the primary component q-^^ 
(such as degree, multiplicity and affine Hilbert function if the component is prime) are the same for 
all the conjugate components cr(qj^''), a G GaZ(L,;/Q), since we actually compute them by an initial 
ideal, which is invariant by conjugacy. 



2 Algebraic extensions of Q and modular computations 

We are interested in preserving some properties of an ideal a in ii = L[X] (L is a normal algebraic 
extension of Q) "modulo" a well-chosen prime integer p. First of all we need to establish how we can 
compute an algebraic number modulo a prime p and then we will see that in general the reduction 
modulo p of the coefficients of a polynomial ideal preserves the affine Hilbert function of the ideal 
itself. 

Let a = (/i, . . . , /s) C L[X] be an ideal. L is a normal algebraic extension of Q of degree 
s: L ~ Q(a), where a is an algebraic number such that its minimal polynomial is q{T) G Q[T], 
degg(T) = s and q{T) = J2i=ii'^ ~ ^il*^)) where Oi are the automorphism of L fixing Q, Oiia) = 
Ui are the conjugates of a over Q. 

In the definition of a reduction of o modulo a well-chosen prime p, our aim is preserving some 
features of a. 

We fix a set of generators (/i, . . . , fs) in Q(a)[X]. We multiply each fi with a scalar q such that 
Ci ■ fi ^ Z[a] [X] and the coefficients of Cj • fi have g.c.d. (on the integers) equal to 1: we call such a 
set of generators in Z[a][X] primitive. We keep on writing fi for Ci ■ fi. 

We now consider a prime integer p such that q{T) sphts in Z/pZ = Fp in the following way: 

q{T) = Si{T)-S2{T)modp, deg Si{T) = l,deg S2{T) = s - 1, gcd(5i(T), 52(T)) = 1. (2.1) 

Thanks to Chebotarev's Density Theorem (1221), we know that there are infinite prime integers p 
for which (12.11) holds. 

Let /3p be the only root of Si{T) in Z/pZ, < /3p < p - 1: 5i(/3p) = 0. We then define the 
following map, from the ring Z[a] of Z to the finite field Fp = Z/pZ: 

tpp : Z[a] Fp 

f3p (2.2) 
a € Z I— )• a mod p. 

This definition on the generators obviously extend to a homomorphism of rings, well-defined 
because of the choice of p and consequently of /3p. 

We can then extend this homomorphism to the polynomials: 

il;p:Z[a\[X.]^Rp = ¥p[X\ 
/ = 5^a/X^^/ = ^Vp(a/)X^ 



7 



If we consider a = (/i, . . . C Q(a)[X], we can assume that the chosen generators are 
primitive and are in l\a\ [X]; we define a = (/i, . . . , /s) C 

Remark 2.1. Observe that the definition of a is independent on the chosen set of generators of a: 
if a = (/i, /s) = (/{,... , //) then (/i, /s) = (/{,... , //) as ideals in Rp. 

Example 2.2. Consider the ideal 

a = (3y2 - 2V3ZX,3YX - V3V2Z,2X'^ - V2Y) C C[X,Y,Z]. 

This set of generators has coefficients in the algebraic extension Q{V2 + \/3), which is normal, 
[Q(\/2 + ^/3) : Q] = 4, the minimal polynomial ofV^ + V^ is q{T) = - lOT^ + 1. 
Consider now p = 23; 

q{T) = {T + 21) • (T + 12) • (T + 2) • (T + 11) mod p. 

We consider the homomorphism ipp such that ipp{-s/2 + -v/3) = 21. With this definition ofipp, we have 
that: ipp{V2V3) = ll,V'p(\/2) = 5,^p(\/3) = 16. 
5o a = (3 y2 + 14 ZX, 3YX + 12Z,2X'^ + 18 Y). 

We can choose a prime p satisfying (12.11 ) using the following lemma. 

Lemma 2.3. [[3], Lemma 12] Let q{T) G Z[r] be a polynomial and p a prime number such that p 
divides q{^), p does not divide the discriminant of q{T) andp > deg(g(T)). 
Then there exists a root in Qp of q{T), considered as a polynomial in Qp[r]. 

We would like to understand in which cases computations modulo a prime integer p preserve the 
affine Hilbert function of the ideal o. From Proposition 11.71 we can bring back our problem about the 
choice of a prime p preserving the affine Hilbert function to the choice of a "good" p preserving the 
initial ideal (w.r.t. some term ordering) of a. 

Assume that (gi, . . . ,gs) is a Groebner Basis of a w.r.t. <, degree compatible term ordering, 
and that these polynomials are primitive; we now choose a prime p satisfying (12.11 ) for the minimal 
polynomial of a. Then, we define a := (51, ... , Qs) ^ Rp. 

Proposition 11.71 gives a necessary condition for a prime integer p to preserve the affine Hilbert 
function of 0: if computations modulo p preserve HF^^^ then they also preserve LM-^(a) with 
respect to a degree-compatible term ordering ^. 

We will show that a finite number of primes p does not satisfy this necessary condition. 

Lemma 2.4. Consider a algebraic number on Q, [Q(a) : Q] = s, Li, . . . , non-zero elements of 
Z[a], Li = X]j=o '^j*'*'^"'- There is a finite number of prime integers p such that ipp(Li) = Ofor some 
i. 

Proof. We will proceed by contradiction. 

Suppose that there are infinite prime integers p such that ipp maps to zero at least one of the Lj's. 
In particular there is an index i such that L,- is mapped to zero by infinite maps il)p. We can define the 
polynomial L{T) = ^^""J afxK 

If deg L{T) = 0, L{T) = Lj is an integer and there is only a finite number of ipp mapping Li to 
zero, because there is only a finite number of p's dividing it. 

We can then assume that d = deg L{T) > 1. We consider a prime p such that 

ML^) = 0, p>\\L{T)rMT)\f, 
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where q{T) is the minimal polynomial of a; we can choose such a p since the set of p's we are looking 
at is supposed to be infinite. 

Observe that ipp{Li) = L{f3p) mod p = 0. This means that both L{T) and q{T) can be divided by 
(T - /3p) modulo p. But since we chose p > ||L(T)p||g(T)||'^, we can apply 121, Lemma 16.20: 

deg(gcd(L(r),(?(r))) >1. 

Since deg L{T) < deg q{T), this contradicts the fact that q{T) is irreducible. □ 

Lemma 2.5. Let a be an ideal in Q(a)[X], ^ a term ordering and G = {gi^ . . . ,17^} a Groebner 
Basis of a with respect to ^, p a prime integer satisfying d2.iD for the minimal polynomial of a. 
Iftpp does not map to any of the coefficients of the leading monomials of the polynomials in G, then 
G = {gi, . . . ,gr} is a Groebner Basis of a with respect to 

Proof. Since G is a system of generators for the ideal 0, then (5 is a system of generators for 0. 
Consider g £ a and its representation with respect to the basis G (eventually multiplying for c G Z, 
to eliminate the denominators): 

r 

g = '^aigi, Oj G R. 

i=l 

Then we have the corresponding representation of g with respect to the basis G: 

r 

g = ^ai5i. 

i=l 

We have that LM{g) = maxj=i^...^r{-^-^^ (oj^j)}, where the "max" is taken with respect to 

Since LM{aigi) = LM{ai)LM{gi) for every i and i/'p does not map to zero the coefficients of 
the leading monomials of the polynomials in the Groebner basis, we immediately have LM [g) G 
{LM{gi), . . . , LM{gs)) and so G is a Groebner Basis for a. □ 

Theorem 2.6. Let a be an ideal in R = Q(a)[X]. Then for a finite number of prime integers p, we 
have that 

HF^^/,{i)^HF-^/,{i) for some i. 

Proof. We fix a degree compatible term ordering ^ and consider the Groebner Basis G = {gi, . . . , , gr} 
of a. Thanks to Proposition 11.71 HF^^^(i) 7^ HF^ ^-(i) for some i only if the initial ideal of a dif- 
fers from 0. 

We apply Lemma 1231 there is only a finite number of primes p such that the initial ideals of a and a 
are different. □ 

Corollary 2.7. There is a finite number of prime integers p such that the affine Hilbert function 
and dimension of and a differ If a is a primary component of some ideal, the same holds for the 
multiplicity of a and a. 

3 An exact strategy: Elimination of Variables and Colon Ideals 

In this section, we will present an exact technique, which combines elimination of variables, absolute 
factorization and computation of colon ideals and gives the primary decomposition of an ideal. 

This technique is intuitive and immediate from the geometric point of view: the elimination of 
variables geometrically corresponds to projection on some linear space. 
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The use of projections reduces the problem of decomposing an algebraic set to a multivariate factor- 
ization; this was first showed at the beginning of the XX-th century (see [16J). After that, Seidenberg 
in EOl used a more rigorous formalism then Hermann's to establish which ideal operation can be 
actually computed, depending also on the features of the polynomial ring we are working on. 
The computation of a colon ideal geometrically corresponds to take off the points of an algebraic set 
from another one. This relation between the "difference" of varieties and colon ideals is well-known 
too (see Q, Chapter 4, §4). 

The strategy of elimination of variables and computation of quotient ideals is also used in ||4|. The 
authors in f4l can move around the computational effort of using generic projections (in the sense of 
Definition 13.21 ) by using another powerful tool, which is the relation between flatness and variation 
of staircases. Thanks to the study of the flatness of the variation of the staircases, they can split the 
ideal according to the splitting of the projection, even if it is not a generic one; they repeat the process 
for the ideals obtained by splitting and at each step they have the dimension or the multiplicity of the 
ideals decreasing, so their algorithm terminates. The algorithm in |4| has different variants and can 
give the strict, isoradical and reduced equidimensional decomposition of a polynomial ideal in Q[X] 
(see m. Introduction, for the different kind of decompositions). 

Unfortunately, even if in f4l the proofs are given for a field k, with suitable properties that Q has, the 
implementation of the algorithm in GoGcAdeals only with Q or a finite field as coefficient fields, so 
this algorithm is not used in a CAS for the computation of an absolute primary decomposition (see 
Section|5]). 

Although the basic idea of the technique we are going to use is well-known, for lack of a complete 
and accessible reference on the whole strategy relating the algebraic and geometrical point of view , 
we will present it in an exhaustive way. 

The decomposition algorithm that we can obtain using these techniques (Algorithm [D is exact 
but not useful in practice: the computations needed are quite long and hard to perform. Anyway, we 
will investigate this method in details since later (Section |4ll we will modify this strategy giving up 
the exact computations in order to gain velocity in computations, but preserving some information 
about the irreducible and reduced components, namely their affine Hilbert Function. 

First of all, we now investigate some properties of the projection of varieties. For the omitted 
proofs, see |[2T1 . 

Definition 3.1. A linear projection is a surjective affine map: 

ft 

P = {xi,...,Xn)eC^^{Li{P),...,Lm{P)) withLi{X) = a,,+Y,ai^Xj. ^^"^^ 

i=i 

In a similar way, we can define a linear projection to projective space, considering 
TTL : \ L ^ 

P=[xo:xi:...:x„]^[Lo(P):Li(P) 

L C P" is the point of intersection of the linear equations Li and it is the center of the projection. 
Two projections vri and 112 from P" to P™ are equivalent if they have the same center L. 

We now consider a projection vr : C" — )• C™ and we restrict it to an equidimensional algebraic 
set W of C" with dimension c. It is not always true that n-p- is proper: for instance, the projection of 



L„(P)], L,(X) 



n 

E 

7=0 



(3.2) 
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the hyperbola defined by the ideal {XY — 1) on the X-axis is the y-axis without the origin; with this 
particular choice of the linear space to project on, the restriction of vr to the hyperbola is not proper. 

Proposition 3.2 (f2T'|, Lemma 5.1). Let W be a closed algebraic set ofC^ all of whose irreducible 
components are of dimension c. For a general linear projection n : C" — )> C"^, m > c + 1, the map 
TTp^ is proper and generically one-to-one. 

In particular ^{W) is a closed algebraic set ofC"^ of degree equal to the degree ofW. 

Corollary 3.3. Let W be an equidimensional algebraic set of dimension c in C". For a general 
projection n : C" — )• C'^'''^, the following hold: 

L ifWi and are two distinct irreducible components ofW, then 'n{Wi) ^ 'n{Wj); 

2. W is irreducible if and only if the polynomial defining tt{W) is absolutely irreducible; 

3. if D{Ti, . . . , Tc) is the polynomial defining the projection Tr{W) and we consider its absolute 
factorization D = ■ ■ ■ -D™', then s is exactly the number of distinct irreducible compo- 
nents ofW, rrii is the multiplicity of the component Wi and deg Di its degree. 

We now assume that W is an equidimensional algebraic set of dimension c defined by the ideal 
C C[X], = {Fi, . . . ,Fr),Fi G Q[X]; we further assume that we performed a generic Unear 
change of coordinates with integer coefficients; in this way we can consider the projections on linear 
spaces defined by equations of kind Xj = to be generic. 

We consider the projection tti : C" — )• Hi, with Hi the linear space of dimension c + 1 defined 
by the equation Xc+2 = • • • = X„ = 0. We call Zi the projection of the algebraic set W on Hi. This 
is a hypersurface in Hi and its decomposition is equivalent to the absolute factorization of the multi- 
variate polynomial Di{Xi, . . . , Xc+i) G Q[X] defining Zi on the linear space Hi. Furthermore, the 
components of Zi are in one-to-one correspondence with the irreducible components of W. 

If Di = • • • e C[X], is the absolute factorization of Di, each factor Dij defines 

in C" a ruled surface, a cylinder, containing the component Wj. 

Algebraically, thanks to the generic change of coordinates and to Corollary [33J if a = P| q^, then there 
is only one absolute factor Dy in the absolute factorization of Di such that D^j^ G q^ . Furthermore, 

for m < ruj, D"^- ^ q^. 

We can then start considering the ideal a + [D^j'). But V{a + {D'^j')) contains not only the 
component Wj, but also the points of the sets V{D^j' ) nW^for k ^ j. 

Thinking of the corresponding varieties ^(qi) and V{D'^j' ), we can classify the primary compo- 
nents of the ideal o + {D^^ ) as follows: 

• a primary component of dimension c, which is the primary component q^ such that S q^ 
(and in this case we write i = j). Obviously in this case q^ + (D^^ ) = c\j and pj + {Dij) = 

• other primary components of dimension < c when V((\i) n V{D^j^ ) is non-empty and it does 
not contain the component Wj . 

Finally, the ideal a+ (D^^ ) defines the component Wj with some extra components of dimension 
< c. In order to avoid these extra components, we can repeat the same steps with (n — c — 1) more 
generic projections vr^ from C" onto the linear spaces H^, defined by equations of kind Xi = 0, 
thanks to the chosen generic coordinates (see also Section 13.11 for more details). If we consider 
= Zk, the polynomial defining Zk on H^ is again a rational multivariate polynomial Dk in 
c + 1 variables. We compute its absolute factorization and we obtain another cylinder containing Wj , 
defined by the factor D^j. 
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Actually, in the absolute factorizations of the polynomials Dk there may be several factors with 
the same degrees and multiplicity. For instance, this happens when one of the components Wj is a 
non-pure rational one (see Definition I1.8I ). This can give ambiguity in matching the factors of the 
polynomials whose zero set contains the component Wj . 

In order to match the factors defining the cylinders containing the same component, we can look 
at the Hilbert dimension of the ideal a + (^fc=i(^^^ ))■ This dimension is c if and only if all of the 
sets V{D^- ) contain the same irreducible component of W; this follows from Corollary 13.31 1. 

In order to find the correct matchings, we will not compute the Hilbert dimensions of the ideals a+ 
(X^feii {P'^j )) fo'' every possible (n — c)-uple of factors; we will compute the Hilbert dimension only 
for the matchings such that all the factors D^j have the same degree and multiplicity. Furthermore, 
we get an almost certain probabilistic check by considering a generic section with a linear space of 
dimension n — c: that is, we can look for the matchings such that the ideal a + {Yl'k=i(^hj' )) ^^^^ 
(n — c) variables evaluated in some integer values is zero-dimensional and nonempty (see Algorithm 

m- 

Once matched the absolute factors found through the projections, {Dkj}j=i^...,c (after re-indexing 
of the factors), the ideal 

c+l 

a,:=a + (^iD^)) (3.3) 
fc=i 

is "almost" the primary ideal corresponding to the component Wj : there are some embedded compo- 
nents left, which geometrically are in Wj n Wi, for I / j. 

Lemma 3.4. Consider (n — c) projections 1:1 from C" to C^^^, pairwise not equivalent. If there are 
points Pi , P2 G such that 

■Ki{Pl) = TTiiP-i) Vz = 1, . . . ,n - c, 

then Pi = P2. 

Proof We write TTk{P) = {Lf\p), . . . with Lf')(X) = a^^ + J2]=i ^If Remark 

(k) (k) 

that since a projection is a surjective map, L;- ^ , 1 < i < j < c + 1; furthermore, we will be 
interested in the matrix whose lines are the vectors ( a,-!^^ ) for A; = 1, . . . n— c, z = 1, . . . , c+ 1. 

V J j=l,...,n 

It is a matrix with (n — c)(c + 1) lines and n columns; since the chosen projections are pairwise non- 
equivalent, this matrix has maximal rank: (n — c) (c + 1) > n for all < c < n — 1, so the considered 
matrix has rank n. 

Pi and P2 have the same image under iTj if and only if L'f\Pi) = L'f\P2), i = 1, . . . , c + 1. 
We obtain (n — c) • (c + 1) equations of kind 

n 

-^j J - Pl-\^j )j=l,...,n,P2 - [Xj )j=l,...,n- 

The unique solution to this system of equations is the trivial one, since the matrix associated to the 
linear system is exactly (^olj''^ whose rank is n, and so we have that Pi = P2. □ 

Corollary 3.5. Consider W equidimensional algebraic set of dimension c in C" and Wi and W2 
irreducible components ofW. If Pi ^Wi and P2 G are such that 

miPi) = 7ri{P2), 

for (n — c) generic projections tTj : C" — )• C'^'^^, then Pi = P2 is a point in Wi Ci W2. 
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Thanks to Lemma [34l and the one-to-one correspondence between the components of W and the 
components of 7r(#'), for a general projection vr, if we consider the ideal aj as in (I3.3I ). V{p-}) is Wy, 
however the polynomials in a^- vanish on the points of n Wi, I / j, with multiplicity higher than 
mult(#;). 

So, for the moment, we have an ideal such that its zero set contains the irreducible component Wj 
but there are some embedded components; we now show that the points of these embedded compo- 
nents are also contained in the zero set of the singular locus of <Xj (Definition l3.6l ). 
Considering again primary decompositions, for what concerns we have that: 



where hi is a primary ideal of dimension < c. Geometrically, the primary components bj correspond 
to the irreducible component of "Wj n j 7^ k. These components are in the singular locus of R/a. 
Here we just recall the algebraic definition of singular locus and the useful Jacobian criterion. 

Definition 3.6. Let a be an ideal in R = ]K[X], IK perfect field, = (/i, . . . , fg)- A prime ideal p 
containing a is in the singular locus of R/a if the localization of R/a at p is not a regular local ring. 

With an abuse of notation, we will say "singular locus of 0" for the singular locus of i?/a. 

Proposition 3.7 (| 1 1 1, Corollary 16.20). Let a be an ideal in K[X.], K perfect field, a is equidimen- 
sional with dimension c, 0= (/i,...,/^). Let J be the ideal generated by the {n — c) x (n — c) -minors 
of the Jacobian matrix {dfi/dXj). Then J defines the singular locus of a: a prime p contains J if 
and only ifp is in the singular locus of a. 

We can then compute easily the equations defining the singular locus of 0: 

• Compute the jacobian matrix of a; 

• Compute {Ml}, the (n — c) x (n — c) minors of the jacobian matrix; 

• The singular locus of is defined by (Mi). 

We are then interested in considering aj and removing the "embedded" primary components 
which contains the singular locus of a. We can do this for the ideals aj containing a reduced compo- 
nent, that is the ones obtained from irreducible factors of multiplicity 1 . 

In general, for ai, a2 ideals in a ring R: 



If we consider an ideal generated by a single element / of R, we will simply write (a : /) instead of 



c+l 



V 




(3.4) 



(ai : 02) = {/ G R\f • aa C ai} 



(a : (/)). 



Lemma 3.8 (H], Lemma 4.4.). Let qbe a p-primary ideal, f an element of R. Then 



L iff G q then (q : /) = (1); 



2- if f ^ <\ then (q : /) is p-primary, and therefore y^(q : /) = p; 



3. iff (/ p then (q : /) = q. 
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Proposition 3.9. In the previous setting, if Dij are factors of multiplicity one in the factorizations of 
the polynomials Di, such Dij G c\j, aj = 0+ (^^jl^(Z)jfcj)) and M is a (n — c) x (n — c) minor of the 
Jacobian matrix of a, then [aj : M) is exactly pj, the prime ideal defining the irreducible component 
Wj ofW = V{a). 

Proof Consider the primary decomposition of aj : 

C+l V 

aJ = a + {Y,{DkJ)) = ^Jr^{[]h,) 

k=l i=l 

where is the only component of dimension c, while hi are embedded components of dimension 
< c. These embedded components correspond to the intersection of Wj with the other irreducible 
components of a. Furthermore, since the factors D^j have multiplicity one in the factorizations of the 
Di, this means (Corollary 13.31 ) that = (\j, that is (\j is a prime component, so we can write pj for 

This means that if M is an equation of the singular locus of a, then M S bj for i = 1, . . . , f and 
so (bj : M) = (1). On the other hand, pj is the only primary (and prime) component of aj which 
does not contain the singular locus: so {pj : M) = pj. 

So 

aj:M = ipj : M) n ^f|(bi : M)^ = p^. 

□ 

Remark 3.10. Proposition 15.91 applies only for components of multiplicity 1 (that is, for factors in 
the absolute factorization of multiplicity 1). 

In fact, if we consider an ideal aj = a + (^^=i{D^-)),mj > 2, we have that this ideal contains the 
singular locus of a, so M ^ aj and (using Lemma [Ol ). (oj : M) = (1). 

For components of multiplicity > 1, we may "clean up " at least some of the embedded compo- 
nents of the ideal aj in the following way: 

• we use the described strategy for the components of multiplicity one, obtaining the set of 
prime ideals {pi^, . . . ,pi„} which are the ideals for the irreducible and reduced components 

• since we cannot use the colon ideal with respect to a generator of the singular locus, we then 
compute 

b:= (...((a,:/,J:/,,)...:/,J, 

for fi^ G piy 

In this way, we can "clean " part of the points of the embedded components. For the remaining 
ones, for any component of multiplicity > 1 we should consider the absolute factors Dij of Di not 
contained in this primary component, and compute again some "nested" colon ideals, starting from 
b, with respect to the Dij. 

3.1 Algorithms 

In this section we will summarize the strategy of elimination of variables and colon ideals, writing 
down the main algorithm for primary decomposition (Algorithm [T|l and an auxiliary one to match 
(n - c)-uples D"^'..., D'^P (Algorithm |2l). 
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In order to consider (n — c) projections vrj : — )• C^^^ (pairwise non-equivalent) and compute 
the polynomials Di which is zero on ■Ki{W), at the beginning of Algorithm [T] we perform a generic 
change of affine coordinates and then use projections on "coordinate" affine linear spaces, namely 
linear spaces defined by equations of kind Xj = 0. While in the previous section we used "H" for 
the coordinate affine linear space, we will now use to denote the set of indexes of the variables 
Xj to eliminate to obtain the projection of the curve: 

Xj=0, jeH = {ji,...,jn-c-i]^l:---:n=[n]. (3.5) 

We can always find (n — c) sets of indexes Hi,i = l,...,n — c, such that Hi ^ Hm for every I ^ m. 
If -H" = {ii, • • ■,in-c-i], we write for {X;^ , . . .,Xi^^^], with [n]\H = {h, . . . ,/c+i}- 

The actual computation of the polynomial Di is obtained by the computation of a Groebner Basis 
with respect to an elimination term ordering. 

If we consider the set of indexes H as in (13.51) . we can compute a Groebner Basis which elimi- 
nates from a the variables Xj, j G H; from now on we will denote such a Groebner Basis with 
G'-Be«im(H)(fl)> or> if we fix Hi, i = 1, . . . , n — c, we simply write GSg;j^(j)(a). Then Di is the only 
generator of Gi?eJim(j) {^) such that 

A e GB^iim(r) (a) n Q[Xi]ie[„]\H, (3.6) 
We will always assume that Hi = {c + 2, . . . ,n}. 

In Algorithm |2l in order to match the absolute factors corresponding to the same component, we 
use the fact that if an algebraic set has dimension < c, then its general section with a linear space of 
dimension (n — c) is empty; this means that if b is an ideal of dimension < c, then the affine Hilbert 
Dimension of b with (n — c) variables evaluated in integer values is generically —1. Furthermore, 
thanks to the generic change of coordinates, we can evaluate the variables in 0. 

Thanks to the results of Section[Tl we can modify Algorithm [T] and [2] and avoid repeating the com- 
putations for conjugate primary components. We present the Algorithms without this for simplicity 
and also because they will not actually be used in this form. 

Indeed, Algorithm [T] is exact, but in practice the computations are too heavy: the main computa- 
tional difficulties are at Step |2] (because of the computation of the Groebner Basis), Step [3] (because 
of the multivariate absolute factorization) and Step [8] (because of the computation of the colon ideal) 
of Algorithm [U 

In the next section we show that we can perform all these computations modulo a well-chosen 
prime integer p, obtaining Algorithm [3] Its output is no more exact, in the sense that it does not 
return the ideals of the reduced components, but it returns information about the ideals of the reduced 
components, in particular their initial ideal and affine Hilbert functions. 

This will be possible observing that the computations of the troublesome steps of Algorithm 
[T] are actually obtained through Groebner Basis and that a good choice of p preserves along the 
computations the initial monomials of the polynomials in the Groebner Basis and so the Hilbert 
functions of the ideals. 
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Algorithm 1: Exact Decomposition of an equidimensional algebraic set 
Data: a € Q[^i, • • • , Xn], a of pure dimension c 

Result: The number of irreducible components, their degrees and multiplicities. 
If the multiplicity of the component is 1, a system of generators of the ideal /(#^); if the 
multiplicity of the component is > 2, the equations of (n — c) hypersurfaces "isolating" the 
component from the other ones. 

1 Preprocessing: Perform a generic affine change of coordinates on the generators of a; 

2 Compute (n — c) polynomials, eliminating from a the set of variables corresponding to the set 
of indexes Hi, i = 1, . . . ,n — c: Di is the only generator of Gi?e/im(i) {^) with the property 
(13.61 ). I = 1, . . . , n — c ; 

3 Perform the absolute factorization of Di, i = 1, . . . , n — c 

4 Match the Dij's through Algorithm |2] in such a way that (after re-numbering of the factors) 
ttj = a + (X]fe=i(^fcj)) contains the component #5-; 

5 Compute the Jacobian matrix of a and a minor of size (n — c) x (n — c), M ; 

6 for ifrom ItosAo 

7 if mj is 1 then 

8 I Compute a^*^^ the quotient ideal of with M: Uj^^^ '■= '■ {M); 

9 end 

10 end 

11 Return: 

12 for j = 1, . . . , s: deg Dij degree of the component, mj multiplicity of the component; 

13 if rrij > 2, aj ideal isolating the component; 

14 ifm, = l, oSf^=/(^). 
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Algorithm 2: Matching of factors through Hilbert Dimension 
Data: a and the absolute factors {Dij}j=i,...,n-c, obtained at Step[3]of Algorithm[T] 

j=l,...,s 

Result: After relabeling the polynomials for the index j, L := [{D-^ )i=i,...,n-c)]j=i,...,s with 
V{<^+ (Z]r=i(-^«i))) containing the component Wj of y(a) for every j = I, . . . ,s. 

1 L := empty list; 

2 for j from 1 to s do 

3 Consider a (n — c)-uple {Dij)i=i^,,,^n-c s-t. deg Dij = dj and niij = mjii; 

4 Compute the Hilbert Dimension h of the ideal a + (X]r=i(-^«j)) '^^^^ {"^ ~ c) variables 
evaluated in 0; 

5 iih=0 then 

6 I add the (n — c)-uple {D'^-^ ) to the list L; 

7 else 

8 I go back to Step[3]and change {n — c)-uple; 

9 end 

10 end 

11 Re-number the factors of Di in such a way that the (n — c)-uples in L are of the form 

12 Return: L 



4 Modular Algorithms 

We use the results of Section |2] on the exact decomposition strategy presented in Section [3] We 
develop an algorithm which takes as input an ideal a with generators in Q[X], defining an equidimen- 
sional algebraic set W in C" of dimension c, and gives as output the number of primary components, 
their degrees, multiplicities and the affine Hilbert function of the components of multiplicity 1 . 

Remark 4.1. We did not present a non-modular version of Algorithm^ since all the algorithms 
presented in Section \3.1\ are not actually used. We insert this further procedure to avoid useless 
computations in the calling of Algorithm\5\in Step^of Algorithm^ 
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Algorithm 3: Modular Algorithm for affine Hilbert Function 



Data: a = (Fi, . . . , Fm), Fi G Q[Xi, . . . , o equidimensional with dimension c 
Result: The degree and multiplicity of the primary components of a; for primary components 
of multiplicity 1, their initial ideal with respect to a degree compatible term ordering. 

1 Preprocessing: Perform a generic integer change of coordinates on a, with coefficients in Z; 

2 Fix (n — c) different coordinate linear spaces of dimension c + 1, defined by set of indexes Hi, 
i = 1, . . . ,n — c; 

3 Compute Di{0,. . . , 0, Xc+i) as the only generator of G'-Beiim(i)(fl|xi=o,...,Xc=o) respecting 
(1X61) : 

4 Compute the rational factorization: Di{0,..., 0, X^+i) = d^^\Xc+i)"'^ ■ ■ ■ d'^P (Xc+i)""'; 

5 for j from 1 to s do 



6 
7 

8 
9 

10 

11 
12 



13 
14 



Choose a prime integer pj dividing dj^^ (0); 

Compute Di mod pj as the only polynomial in GBgHfy^^^i^ (a) mod p respecting (13.61 ). 
i = 1, . . . ,n — c; 

Compute the modular factorizations Di mod pj, i = I, . . . ,n — c; 

Apply Algorithm m to the rational factor d^^^ [Xc+i) and the modular factors of Di, 

obtaining the set of modular factors Vj of Di ; 

Choose D\ mod p, G Vj of minimal degree such that r,- = ^ G 

jk J J degO'-V 



jk 



if Vj > 2 then 



Apply Algorithm [5] to match the modular factor D^^~^ with D^^l modular factor of Di, 

jk jk 

obtaining a^-* = a + Y17=i (-^^2)"^^ mod pj (after re-labehng of the factors); 

JK 

if rrij is 1 then 

Compute the Jacobian matrix of a mod pj and a minor of size (n — c) x (n — c), 
Al- 
ls Compute (aj''^ : M) mod pj; 

16 end 

17 end 

18 end 

19 Return: s number of rational components 

20 for every j = 1, . . . , s 

21 Vj number of non-rational components constituting the rational component qj 

22 deg D^^^ degree of each non-rational component of qw 

jk 

23 TTij multiplicity of the non-rational component 

24 if rrij = 1, pj and (af^ : M) mod pj ideal having the same initial ideal and same Hilbert 
function as qij 

25 if rrij > 2, p and D^^^, i = 1, . . . ,n — c, image modulo p of (n — c) a polynomial contained in 

jk 

ql but not in the other absolute components of a; 
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Algorithm 4: Partition of modular factors 



Data: d{Xc+i) G Q[Xc+i], an integer p dividing d{0) and 

D{Xi, Xc+i) = ULi A(Xi, . . . , Xc+ir^ mod p such that 
d{Xc+i) I D{0,...,0,Xc+i) modp 

Result: A set containing the modular factors of d{Xc+i) mod p 

1 A := empty list, i := 1, 6 := 1; 

2 while i < / do 



3 

4 

5 

6 

7 

8 

9 
10 
11 

12 end 

13 Return: A 



nii := multiplicity of d{Xc+i) in the rational factorization of D{0, • • • , 0, Xc+i) 
if Di{0, . . . , 0, Xc+i) mod p divides d{Xc+i) mod p then 



add Di{Xi,... 

if<5 = d(Xe+l) 

I i:=l + l; 
end 

end 



,Xc-\-i) mod p to ^; 
.,0,Xc+i) modp; 
mod p then 



Algorithm 5: Matching of modular factors through affine Hilbert Dimension 
Data: D^^^ modular factor of Di{Xi, . . . , Xc-i-i) and {D^f^^ mod p}j=2,...,n~c modular 

k=l,...,m 

factors of Dj 

Result: aj = a + (-Dp^)™' mod p with Hilbert dimension c 

1 Consider a (n — c — l)-uple (-D["'^)j=2,...,n-c such that deg(-D|,''^) = deg(Dj^^'') and = m, ; 

2 Compute h =Hilbert Dimension of a + + X]"=2 (-^If^)™'' mod p with {n — c) variables 
evaluated to 0; 

3 if /i = then 

I (i) (i) 

4 renumber the modular factor putting := 

5 else 

6 I go back to Step[T]and change (n — c)-uple. 

7 end 

8 Return: o^' = o + {J2]Zt D-^^)""' mod p. 
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4.1 Proof of Algorithm |3] 

We apply the results of Section |2] to the decomposition strategy explained in Section [3] and to the Al- 
gorithms of Section im Again, we deal with an equidimensional polynomial ideal a = (Fi, . . . , Fm) 
with dimension c, Fj G Q[X]. The key point of Algorithm [3]is the choice of a prime integer pi which 
gives a "modular image" of the algebraic number s.t. q.''^ C Q(aj)[X]. For all the notations used, 
we refer to Section [TTl Algorithms [3j HI and [5j 

We will now follow the steps of Algorithm |3] in order to show that it gives a correct output. 

In Step [T] of Algorithm m as in Algorithm [T] we perform a generic change of coordinates; thanks 
to this, the projections on the "coordinate" linear- spaces of dimension c + 1 are "generic" in the sense 
of Proposition 13.21 the components of the projected algebraic set are in one-to-one correspondence 
with the components of the algebraic set itself (see also Corollary 13.31 ). Furthermore, consider on 
one hand the absolute factors of the polynomial whose zero set is the projected algebraic set and 
on the other one the primary components of the ideal defining the algebraic set: factors and primary 
components are in one-to-one correspondence and the degree and multiplicity of a factor is the degree 
and multiplicity of the corresponding primary component (in the sense of Definition 11.51 and 11.61 ). 

In Step[2]we fix (n — c) distinct "coordinate" linear spaces Hi (as explained in Section [37T] ). Using 
projections on these linear spaces, we would like to apply the techniques for absolute factorization 
developed in [3j, but we have to be careful because we do not have one of the main hypothesis: the 
Input of the Abs-Fact Algorithm presented in 131 is a rationally irreducible polynomial. This is not 
our case, this is why in Step 3 of the algorithm we compute a univariate factorization. 
Indeed, assume that we are able to compute Di, the only polynomial in the first c + 1 variables of 
G'-Be««m(i)(ci)- This multivariate polynomial in general is not rationally irreducible; furthermore it is 
not advantageous to compute the multivariate rational factorization of Di. 

We rely on Hilbert's Irreducibility Theorem: for infinite integer specialization of c variables, a ratio- 
nally irreducible factor of the polynomial Di stays rationally irreducible. This means that if 

Z)i(Xi, . . . = . . . • • • . . . G Q[Xi, . . . 

then for infinite xi, . . . , G Z the rational factorization of . . . , Xc, Xc+i) is exactly 

. . . , Xe, Xe+l) = d!i\xi, ...,Xc, X.+i)^' ■ ■ ■ ...,Xc, X.+i)'"^ G Q[X,+i]. 

Thanks to the generic change of coordinates of the Preprocessing Step, we can take xi = • • • = 
Xc = 0. In order to compute this rational univariate factorization without computing Di {Xi , . . . , Xc+i ) , 
in Step |7] we simply specialize c variables of Fi, . . . , Fm and then compute the elimination Groebner 
Basis: 

Di{0,. . . ,0,Xc+i) G GBeiim{i){a\xi=o,...,Xc=o)- 

Since we are considering a generic projection, a rational factor of Di{Xi, . . . , Xc+i) corresponds 
to a rational component of the algebraic set W = V{a) (in the sense of Definition 1 1.8 1 ), while each 
absolute factor of corresponds to an irreducible component. 

Once computed in Step|4]the univariate rational factorization, we then proceed in order to "break" 
the non-rational components. 

We consider the j-th factor of the rational factorization of Di, that is d^^\Xc+i) which has multi- 
plicity rrij. If the corresponding factor d^p {Xi, . . . , Xc+i) of the univariate rational factorization of 
Di{Xi, . . . ,Xc+i) is not absolutely irreducible, then its absolute factors have coefficients in some 
algebraic extension Q{aj). Using [3|, Lemma 11, we can assume that the algebraic extension Q(aj) 
is generically generated by the evaluation of one absolute factor in a point with integer coordinates. 
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Thanks to the generic change of coordinates, we will choose (0, . . . , 0) G 

We choose an integer prime pj dividing dj (0) (Step ^ applying Lemma 12.31 and, relying on ran- 
domness, we assume that the chosen prime pj will preserve the initial ideal of the Groebner Basis 
we will compute along the "FOR" loop (as in Lemma 12.51) . Thanks to Lemma 12.31 if we factor 
Di {Xi , . . . , Xc-\-i ) modulo this prime pj , the rationally irreducible factor dj^^ {Xi , . . . , Xc+i) splits 
(if it is not absolutely irreducible). The homomorphism ipp^ of (12.21 ) is implicitly defined. 
Actually we do not compute Di{Xi, . . . , Xc+i): in fact, in Step |7] we compute directly the modular 
elimination Groebner Basis and then the modular factorizations (Step [8]). In Step |9] we group the 
modular factors corresponding to d^p (Xc+i) using Algorithm^ 

If the rational factor d^^\Xi, . . . ,Xc+i) is absolutely irreducible, then it does not further split 

modulo Pj, that is rj, the number of modular factors of d^^ , is exactly 1. In this case, we can stop 
here and repeat the loop for the next rational factor. 

If dj^\Xi, . . . ,Xc+i) is absolutely reducible, then rj > 2 (thanks to the choice of pj according 
to Lemma [23] ): in Step [10] we choose a modular factor among them having minimal degree which 
divides deg d^-^^ (Xc+i); we assume that this factor is D^^^{Xi, . . . , Xc+i). 

In Step [12] we look for the corresponding modular factor of Di, i = 2, . . . , re — c. Using Algorithm 
[5] we obtain the ideal a'- = a + Y^i=i {D^^)"^^ mod pj with Hilbert dimension c. Corollary 12.71 

k jk 

certifies that o- = thn ia'-)- 

Once defined in Step [12] the ideal o| (re-ordering the indexes) with affine Hilbert dimension c, if 
the multiplicity rrij is 1, we can keep on following Steps [5] and [8] of Algorithm [T] we compute the 
Jacobian Matrix of a mod pj and consider one of its (re — c) x (re — c)-minors, M. We compute the 
colon ideal of a| with M. Let M be the (re — c) x (re — c) minor of the Jacobian matrix of a s.t. M 

mod Pj = M. 

We need to show that for infinite primes pj the colon ideal modulo pj has the same affine Hilbert 
function of the colon ideal in Q(aj)[Xi, that is ■0^^ (a^ : M) = (o^ : M). 

First of all, observe that ai and the corresponding non-modular ide al ai have the same Hilbert func- 
tion for all but a finite number of prime integers (thanks to Theorem l2.6l) . 

Furthermore, we can assume that we compute Jacobian matrix of o and a minor M and then reduce 
modulo Pj. For what concerns the colon ideal, the actual computation is performed using a Groebner 
Basis (see [T|, Chapter 4, §4, Theorem 1 1 for the details). This means that again we can apply Lemma 
I2.5l and so there is only a finite number of primes pj such that V^p^ (o| : M) and (a| : M) differ. 

Remark 4.2. Actually, Algorithm\3\is a Las-Vegas one, just like the Abs-Fact Algorithm of [3]: in 
fact, in the Preprocessing Step, we have to assume that the coefficients for the generic change of 
coordinates are taken in a finite set S £ Z. 

We shall then modify the Preprocessing Step of Algorithm\3\and insert a small loop in SteplWl in 
order to stop the execution and go back to the Preprocessing Step, if we cannot define a rj G Z (see 
Abs-Fact Algorithm of [3]). 

We can also compute the "minimal " rational algebraic extension Lj = Q(aj) containing a 
set of generators of the ideal Oj. We can apply the LLL method developed in [3]. Unluckily, we 
do not have a technique to estimate the needed level of accuracy. We can just try to compute the 
minimal polynomial which defines Lj with increasing levels of accuracy and stop when we get the 
same polynomial q(T) with 2 different levels of accuracy. 
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5 Tricks on an example 



The data and Maple files of the examples we are going to discuss are available at 
https://sites.google.com/site/cristinabertone/examples-for-modular-decomposition 

We now test our algorithm on a quite simple example (see the file DecompositionCIcurveDe- 
gree48.mw). 

We consider a complete intersection ideal a C C [X, Y, Z] generated by two polynomials with 
rational coefficients, F,G ^ Q[X, Y, Z], of degree 8 each, rationally irreducible. 
We consti^ucted this c.i. curve in such a way that we know that it has non trivial primary compo- 
nents, in particular it has a rational primary component of degree 14, that splits in 2 absolute primary 
components of degree 8 each, generated by polynomials in Q{y/2)[X, Y, Z\. 

The complete intersection curve C = V{a) has degree 48 (one can see this, for instance, using a 
generic plane section and counting points with multiplicity). 

Since a is generated by 2 polynomials, we can use resultants instead of elimination Groebner Basis 
to compute the elimination of variables. We perform a generic linear change of coordinates and we 
compute 

r :=i?esz(i^(0,y,Z),G(0,y,Z)), 
which has degree 48 and factors over the rationals (in less than 1 second) as: 

• (y) factor of degree 14 and multiplicity 1; 

• dg^-* (y) factor of degree 4 and multiplicity 1; 

• (y) factor of degree 22 and multiplicity 1; 

• ^4^'' (y) factor of degree 2 and multiplicity 2; 

• d^l^ (y) factor of degree 1 and multiplicity 4. 

So, using Definition 1 1.81 the complete intersection a has 5 rational components q^, three of them with 
multiplicity 1, with degrees given by degdj-^^(y) (thanks to Corollarv 13.31) . 

We can proceed in the following way: for each rational factor (y), we choose a prime number 

'Pi dividing (i.^^(O), except for i = 5: indeed, we do not look for a prime dividing ^^^^(O), since this 
rational component will not further split. We then compute the projections on the coordinate plane 
modulo Pi, i = 1, ... ,4. Then we compute the modular polynomial describing the projection of 
the curve for each prime pi and its modular factorization. We know that pi forces the rational factor 
corresponding to d[^^ {Y) to split (if it is absolutely reducible). We check whether there is a prime pj 
between the chosen ones such that it forces all of the rational factors: if we find one, we can perform 
all of the computations modulo this prime. If not, we can in any case choose to compute modulo 
"some" of the primes pi: if pi and pj both give the desired splitting for the i-th and j-th rational 
factor, then we can compute the corresponding ideals modulo pi (and not use pj). 

For the ideal a, we see that p = 89 give the desired splitting for all of the 4 rational factors which 
may be absolute reducible. So we will compute only modulo 89. 

The computations of Di{X, Z) = ResY{F, G) mod p, D2{X, Y) = Resz{F, G) mod p and 
their modular factorization take less than 4 seconds each. 

All of the rational factors further split modulo p in 2 factors. So we compute the initial ideal 
(and affine Hilbert function) for one of the two absolute components of degrees 7, 2 and 11. The 
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other absolute components have multipUcity > 1, so we do not perform on them Steps [T4l and [T5] of 
Algorithm [21 

Performing Algorithm |3] (including the matching of the factors through Algorithm [S]), we then 
obtain the initial ideals: 

. deg q^^^ = 7, mueM?) = {X\Y\ X^Z\ X^Y, XY\ X XY^ XY Z^); 
. degq^^^ =2,intuM^) = {X,Y^); 

• degg^^^ = 11, intiexiqf^) = 

= {X^,Y^^,XZ^, X^Y^,XY^, X^Z^,XY'^Z, X^Y Z"^ , XY^ Z^ , XY^Z^,XYZ'^). 

Finally, we use the techniques of the Abs-Fact Algorithm of [3] to compute the polynomial 
q{T) G Z[T] which defines the algebraic extension containing the coefficients of a set of genera- 
tors for the absolute primary components. However we do not have an a priori bound on the size of 
the coefficients of q{T), as pointed out in Remark l42l ). We perform an Hensel Lifting of a modular 
univariate factor until a quite high level of accuracy (in this case, until p^^^); we then construct dif- 
ferent candidates for the minimal polynomial, starting with accuracy p^^, until two different levels of 
accuracy give the same polynomial. 

For accuracy p^^, we see that the minimal polynomial "stabilizes": 

qiT) = 26301054375 - 214355874045600 T + +436754388124393216 

Obviously, since degq{T) = 2, we can easily find a better presentation of the extension Q(q;) com- 
puting the roots of q{T): we obtain that the extension of Q we need can be generated by \/2. 

Summing up, we obtained that the complete intersection curve o = (F, G) C Q[X, Y, Z] has the 
rational primary decomposition 

= qi n q2 n qs n q4 n q5 

with degqi = 14, degq2 = 4, degqs = 22, degq4 = 2 and degq5 = 1 and multiplicities mi = 
iT^2 = ''7^3 = 1, ni/^ = 2, ms = 4. 

Each of the rational primary ideals with multiplicity 1 further decomposes as 



with q^^ C Q{V2)[X,Y,Z], qf^ = a (qf^), where a(^/2) = -^2. 

The whole computation took less than 15 minutes on a home-use personal computer, without any 
problem with memory allocation. 

We point out that it is not that obvious to obtain this kind of information about the decomposition 
of the ideal a. For instance, one may use one of the most popular Computer Algebra System, Maple 
[181. 

We tried to use the Maple command PrimaryDecomposition (whose algorithm is based on 
lfT3l ). which gives as an output the primary decomposition of the ideal a. As input, we also gave 
the algebraic extension of the rationals in which one can find the generators of the absolute pri- 
mary decomposition of a, namely Q{\/2). Even with this further information about the decom- 
position (which is not a priori known from the only knowledge of the rational generators of a), 
PrimaryDecomposition in Maple caused a problem with memory allocation (reaching about 
2.3 GB), after computing for more than 1 hour. 

For what concerns other computer algebra systems, we also tried Singular ([14]), another computer 
algebra system for polynomial computations. We tried to obtain the rational primary decomposition 
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of a using primdecGTZ and the primary decomposition over Q[X, Y, Z] using absprimdecGTZ 
(which are based on [13], the algorithms are described in [9|). In both cases we stopped the compu- 
tations after 2 hours, without obtaining the primary decomposition. 

The CAS CoCoAdSl) has a command PrimaryDecomposition to decompose only monomial 
square-free ideals. It also has a command called EquilsoDec which computes an equidimensional 
isoradical decomposition of o, i.e. a list of equidimensional ideals bi, . . . , bfe such that the radical of 
a is the intersection of the radicals of bi, . . . , bfc. This command is based on the algorithm presented 
in m and it works only using Q or finite fields as coefficient ring. Nevertheless, Equi I soDec could 
not give the output on our example after more than two hours computing. 

Although for the moment we cannot really compare our algorithm with the above ones mentioned, 
we can see that the problem we are facing is challenging and that our modular strategy may move 
around the computational problems of primary decomposition. Nevertheless, we cannot give com- 
plete comparaison for the moment, since we cannot compute the complete primary decomposition 
nor a reduced decomposition: using Algorithm [3] we get several interesting data about the absolute 
primary components of an equidimensional ideal. Indeed, these data may be useful as a guide or a 
bound for numerical algorithms, such as the ones in |[T2ll or ll2T1l . 

Furthermore, the technique taken from Q to construct the algebraic extension containing the co- 
efficients of a set of generators of a primary component, could be used with other algorithms: for 
instance, one can see the example DecompositionCIcurveDegree36.mw . With our modular strategy, 
we obtain not only the initial ideal of all of the absolute primary components of a complete intersec- 
tion in Y, Z] generated by 2 polynomials both of degree 6, but also that Q(\/2) contains the 
coefficients of the generators of these components. 

If we ask Maple to compute the primary decomposition of this ideal, again it does not reach the result 
after one hour of computation. But if we pass to the command PrimaryDecomposition also 
the information that the primary components are inside <[^{^/2)[X, Y, Z], then we obtain the primary 
decomposition in less than 3 minutes. 

This example suggests that we may combine at least a part of our strategy (for instance, the 
construction of the "splitting" field) with an existing algorithm (such as the one implemented in 
Maple), in order to obtain a complete result. 

Conclusions and future work 

In this paper we designed an algorithm which, given a set of polynomials with rational coefficients 
defining an equidimensional ideal a, returns the initial ideal of each absolute prime component of 
a. Furthermore it also returns information concerning non-reduced primary components, such as 
their number, degree, multiplicity. The main ingredients of the algorithm are the classical technique 
of projection and the use of computations modulo well-chosen primes, as done in |31 in order to 
decompose a bivariate polynomial. The obtained results seem promising, mostly for what concerns 
complete intersections and more precisely, for curves in C^. 

A further step is to implement this algorithm in a Computer algebra System, for instance in 
Mathemagix [19]. Mathemagix is a free computer algebra system under development, which has 
available libraries for algebraic computation (such as large numbers, polynomials and others) for 
exact and approximate computation. This should make Mathemagix particularly suitable as a bridge 
between symbolic computation and numerical analysis. 

Our final aim is actually to design an algorithm which uses projections with modular techniques 
and can return the complete absolute primary decomposition of the ideal given as input. The main 
obstacle to this is the absence of a tool similar to Hensel's Lifting (see for instance Il23l . Chapter 15, 
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Section 4), which allows lifting a modular factorization to a rational one; we would need a general- 
ization of this in order to lift the modular decomposition of an ideal. 

Our next task is then to develop such a tool, design and implement a primary decomposition 
algorithm and compare its efficiency with other implemented routines. With this tool, we may be able 
to obtain a complete decomposition for the reduced part of an equidimensional ideal. 

Other possible improvements of the algorithm are dealing with a non-equidimensional ideal a and 
computing also non-reduced primary components; for the first part we will have to choose whether 
dealing only with the top-dimensional part of o or studying also smaller components; for the second 
one, if we will be able to compute these non-reduced primary components, then we will also be able 
to compute a radical decomposition (see [4|). 

We aie hopeful that our techniques are competitive, since at this moment we can already get a lot 
of information concerning the absolute decomposition of an ideal, in a reasonable time with a limited 
use of memory, while other CAS cannot really deal with the primary decomposition of the same ideal 
(see Section [5]). So, even if at the moment the results of our method are partial and cannot be directly 
compared to the performances of other softwares, we believe that this method is on the right path to 
get an efficient primary decomposition algorithm. 
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